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ABSTRACT 

Sgr A* is currently being fed by winds from a cluster of gravitationally bound young 
mass-loosing stars. Using observational constraints on the orbits, mass loss rates and 
wind velocities of these stars, we numerically model the distribution of gas in the 
~ 0.1-10" region around Sgr A*. We find that radiative cooling of recently discovered 
slow winds leads to the formation of many cool filaments and blobs, and to a thin and 
rather light accretion disc of about an arcsecond scale. The disc however does not ex- 
tend all the way to our inner boundary. Instead, hot X-ray emitting gas dominates the 
inner arcsecond. In our simulations, cool streams of gas frequently enter this region on 
low angular momentum orbits, and are then disrupted and heated up to the ambient 
hot gas temperature. The accreting gas around Sgr A* is thus two-phase, with a hot 
component, observable at X-ray wavelengths, and a cool component, which may be 
responsible for the majority of time variability of Sgr A* emission on hundred and 
thousand years time-scales. We obtain an accretion rate of a few xl0 -6 M Q year -1 , 
consistent with Chandra estimates, but variable on time-scales even shorter than hun- 
dred years. These results strongly depend on the chosen stellar orbits and wind pa- 
rameters. Further observational input is thus key to a better modelling of Sgr A* wind 
accretion. 

Key words: Galaxy: centre - accretion: accretion discs - galaxies: active - methods: 
numerical - stars: winds, outflows 



1 INTRODUCTION 

Sgr A* is identified with the Mbh ~ 3 x 10 6 M super- 
massive black hole (SMBH) in the centre of our Galaxy (e.g., 
Schodel et al., 2002; Ghez et al., 2003). By virtue of its lo- 
cation, Sgr A* may play a key role in the understanding 
of Active Galactic Nuclei (AGN). Indeed, this is the only 
AGN where recent observations detail the origin of the gas 
in the immediate vicinity of the SMBH capture radius (e.g., 
Najarro et al., 1997; Paumard et al., 2001; Baganoff et al., 
2003; Genzel et al., 2003b). This information, missing for all 
other AGN because of too great a distance to them, their 
large luminosity, or both, is absolutely necessary if the ac- 
cretion problem is to be modelled self-consistently. 

Arguably the most famous puzzle of Sgr A* is its low 
luminosity with respect to estimates of the accretion rate 
at around the capture radius, i.e. at distances of order 
1" ~ 10 5 7?s ~ 0.04 pc, where Rs is the Schwarzschild radius 
of Sgr A* . Two methods have been deployed to obtain these 
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estimates. From Chandra observations of the Galactic centre 
region, one can measure the gas density and temperature 
around the inner arcsecond and then infer an estimate of 
the Bondi accretion rate of M ~ 10 -6 Mq year -1 (Baganoff 
et al., 2003). However, unlike in the classical textbook prob- 
lem (Bondi, 1952), hot gas is continuously created in shocked 
winds expelled by tens of young massive stars near Sgr A*, 
and there is neither a well defined concept of gas density and 
temperature at infinity, nor one for the gas capture radius. 

The other method addresses this problem by direct 
modelling of the gas dynamics of stellar winds, assuming 
that the properties of the wind sources are known. Three 
dimensional simulations of wind accretion around Sgr A* 
were performed by Coker & Melia (1997), who randomly 
positioned ten mass-loosing stars a few arc-seconds away 
from Sgr A*. They presented two different runs in which 
the stars were distributed in either a spherically isotropic or 
a flattened system. Rockefeller et al. (2004) used a particle- 
based code with more detailed information on stellar co- 
ordinates and wind properties. However, in both cases the 
stars were at fixed locations, whereas in reality they fol- 
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low Keplerian orbits around the SMBH. The accretion rate 
on to Sgr A* predicted by both studies was estimated at 
~ few x Kr 4 M year -1 . Finally, Quataert (2004) stud- 
ied the problem in the approximation that there is an infi- 
nite number of wind sources distributed isotropically around 
Sgr A* in a range of radii. His model yields an accretion 
rate estimate of ~ few x 10~ 5 M Q year -1 . 

Due to recent impressive progress in the observations 
of Sgr A*, we now know much more about the origin of 
the gas feeding the SMBH. Stellar wind sources are locked 
into two rings that are roughly perpendicular to each other 
(Paumard et al., 2001; Genzel et al., 2003b). In addition, the 
wind velocities of several important close stars were revised 
downward from ~ 600 km s -1 (Najarro et al., 1997) to only 
~ 200 km s -1 (Paumard et al., 2001), making the Keplerian 
orbital motion much more important. 

Motivated by these points, we performed numerical sim- 
ulations of wind accretion on to Sgr A* including optically 
thin radiative cooling and allowing the wind-producing stars 
to be on Keplerian orbits. In this Letter we report our most 
important findings: (i) the slow stellar winds are susceptible 
to radiative cooling once they are shocked in the hot bubble 
inflated by the fast winds; (ii) a disc/ring is formed on ~ 
arcsecond scale; (iii) cool gas blobs frequently enter the inner 
arcsecond on low angular momentum orbits, are torn apart 
and thermalized, and then mixed with the hot gas there. Gas 
is thus distinctly two-phase in the inner region, with the cold 
phase being invisible to Chandra, (iv) The resulting 'accre- 
tion rate' (see §4) is of order several xl0~ 6 M Q year -1 , con- 
sistent with Chandra estimates, but is highly variable. This 
warrants a somewhat worrying question: how representative 
is the current low luminosity of Sgr A*? Summarising these 
points, it appears that the accretion flow on to Sgr A* and 
other low luminosity AGN cannot be fully understood based 
on observations of hot X-ray emitting gas alone. 



2 ANALYTICAL ESTIMATES 

The density of hot gas 1.5" away from Sgr A* is about n c = 
130 cm -3 (Baganoff et al., 2003), and the gas temperature is 
T g « 2 keV. The pressure of the hot gas is thus p th — n c T g w 

3 x 10 9 K cm . The temperature resulting from collisions 
of stellar winds with wind velocity « w = 10 s vs cm s _1 is 

T = 1.2 x 10 7 no. 5 vl Kelvin, (1) 

where /10.5 is the mean molecular weight in units of half the 
proton mass. This temperature is to a large degree compat- 
ible with the T g ~ 1.2 keV measured by Chandra at slightly 
larger radii. The optically thin cooling function, dominated 
by metal line emission, is A w 6.0 x 10" 23 Tf °' 7 (Z/3), where 
T7 = T/10 7 K, and Z is the metal abundance relative to 
Solar (Sutherland & Dopita, 1993). The cooling time is thus 

• 3fcT 4 !7 5.4 _! . . 

001 = ~An Mo ' 5 Vs P * y6arS ' ' 

where p, = nT/p t h- The dynamical time near this location 
is i dyn = R/vk — 60 (R") 3 '" 2 years, where R" is the radial 
distance to Sgr A* in arc-seconds and i>k is the Keplerian 
velocity at that distance. This shows that cooling is of no 
importance for the gas originating in the winds with out- 
flow velocity t> 8 ~ 1. However, if the wind velocity is, say, 



t> w = 300 km s~\ then t coo i is only roughly 15 years, which is 
shorter than the dynamical time. Therefore one could expect 
that slower winds may be susceptible to radiative cooling in 
the high pressure environment of Sgr A*. In reality the gas 
velocity is the sum of the wind velocity and the stellar or- 
bital motion. At a distance of 2", for example, the Keplerian 
circular velocity is about 440 km s _1 . With respect to the 
ambient medium, the leading wind hemisphere will move 
with velocity vs ~ 1, whereas the opposite one will move 
even slower than v w . This shows that even for winds with 
velocities « w £ 500 km s _1 , lagging regions of the wind may 
still be affected by cooling. 



3 METHOD AND INITIAL CONDITIONS 

A full account of our numerical method along with vali- 
dation tests will appear in Cuadra et al. (2005). Here we 
only briefly describe the method. We use the SPH/iV-body 
code GADGET- 2 (Springel et al., 2001; Springel & Hern- 
quist, 2002) to simulate the dynamics of stars and gas in 
the gravitational field of the SMBH. This code, developed 
for cosmological simulations, takes into account the (Newto- 
nian) TV-body gravitational interactions of all particles and 
also follows the hydrodynamics of the gas. We use the cool- 
ing function cited in §2 with Z — 3, and set the minimum 
gas temperature to 10 4 K. 

We model the SMBH as a heavy 'sink' particle (Springel 
et al., 2004; Di Matteo et al., 2005), with its mass set to 
3.5 x 1O 6 M . For scales of interest, the black hole gravity 
completely dominates over that of the surrounding stars and 
gas. The inner boundary condition is specified by requiring 
the gas passing within the radius Ri n from the SMBH to dis- 
appear in the black hole. In addition, particles at distances 
larger than the 'outer radius' R out are of little interest for 
our problem and are simply eliminated. 

Stars are modelled as collisionless particles moving in 
the potential of Sgr A*. The stars emit new gas particles 
that are initialised with the minimum temperature and a 
mass mspH = 5 x 10~ 7 M©. The initial particle velocity 
is the sum of the orbital motion of the star, and a random 
component. The latter is equal in magnitude to the wind 
velocity and its direction is chosen randomly to simulate 
isotropic winds in the frame moving with the star. 

Following results of Paumard et al. (2001), we assume 
that 'narrow line stars' produce winds with velocity «„ = 
300 km s _1 , whereas the 'broad line stars' produce winds 
with t> w = 1000 km s _1 . We refer to the former as LBV 
stars, and to the latter as WR stars. The radial extent of the 
inner stellar ring is from 2" to 5", and the inner and outer 
radii of the outer ring are 4" and 8", respectively. The rings 
are perpendicular to each other for simplicity (Genzel et al., 
2003b, conclude that the rings are inclined at i « 74° to each 
other). We argue that the total number of wind sources is 
likely to be higher than those that have been resolved so far. 
On average they would clearly have to be less powerful than 
estimated by Najarro et al. (1997). We therefore use 20 wind 
sources in total, with each star loosing mass at the rate of 
M* = 4 x 1O~ 5 M year^ 1 . Note that this is still a factor 
of ~ 2 — 3 below the observationally estimated total mass 
loss rate from the Sgr A* star cluster, but some of the wind 
sources are likely to be outside R out and thus should not be 
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included here. Typically, the time-step of the calculation is 
~ 0.18 years, so the mass loss rate above implies that ~ 20 
SPH particles are created around each star per time-step. 
The stars are distributed in rings uniformly in radius but 
randomly along the azimuthal angle. Stars in the same ring 
rotate in the same direction, of course, as they follow circular 
Keplerian orbits. We populate the inner ring with 6 LBVs 
and 2 WRs, and the outer one with 3 LBVs and 9 WRs. 

To increase the resolution in the inner region, we split 
the SPH particles that get closer than 10 Ri n to the SMBH. 
To avoid numerical problems, the splitting is performed at a 
randomly chosen time - of the order of the dynamical time 
at that radius - after the particle entered the inner 10-Rm- 
Then the particle is divided into N sp nt = 5 new ones that are 
placed randomly within the smoothing length of the original 
one. The mass of the old particle is equally divided between 
the new ones, while the temperature and velocity are kept 
constant. The rest of the SPH properties are updated self- 
consistently by the code. Other simulations, without split- 
ting, show the same overall results as the one presented here 
(Cuadra et al., 2005). 

The values for the inner and outer boundary conditions 
are i?j n = 0.07" and R ou t = 9". We start with a negligible 
gas density inside the computational domain and then fill 
it up self-consistently with stellar winds. We ran the sim- 
ulation for about 4,000 years. Realistic longer simulations 
would have to include a highly complex gas density field on 
larger scales around Sgr A* and frequent supernovae (see 
below), which we leave for future investigations. The total 
number of SPH particles, JVsph, initially increases steeply 
with time until the hot winds fill up the whole inner sphere. 
Later, the number of particles continues to increase due to 
the fact that most of the slow wind particles are bound. By 
the end of the simulation we reach Asph = 1.7 x 10 6 . 



4 RESULTS 

4.1 Gas morphology. 

Figure 1 shows the column density of gas and the stellar 
wind sources in the inner 6" of the computational domain. 
The inner and outer rings are viewed at inclination angles 
of 40° and 50° in the Figure, respectively. The inner ring ro- 
tates clockwise in this projection, while the outer one rotates 
in the opposite sense. 1 

As fast stellar winds fill the available space, the slower 
ones are shocked and cool radiatively. Dense shells are 
formed around the shocks. These shells are torn into fila- 
ments and blobs. Different parts of the shells have different 
velocities and thus angular momenta, creating many cool gas 
blobs with different velocities. Some are directed to outer 
radii and have velocities large enough to escape from the 
computational domain. Others receive velocities directing 
them inward. At smaller radii the density of blobs is higher, 
thus they collide, and settle in the plane established by the 
inner LBVs. The blobs are also sheared by the differential 
Keplerian rotation, and a gas disc is born at radii somewhat 
smaller than that of the stellar orbits. The blobs typically 

1 An animated movie of this simulation is available at 
http : //www. mpa-garching.mpg.de/ ^jcuadra/Winds 
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Figure 1. View of the gas column density and of the stars (over- 
layed) in the inner 6" cube of the computational domain, after 
about 3,000 years. The 'LBV stars within the slice are marked 
as stars, while the 'WR' stars are shown as diamonds. Stars 
painted in red and black belong to the inner and the outer rings, 
respectively. 



consist of a few hundred SPH particles, many more than the 
minimum 40 neighbours used in the SPH kernel averaging 
for these simulations. 

We find that the disc does not extend all the way in- 
wards, as the inner arcsecond or so is dominated by the hot 
X-ray emitting gas. Surprisingly, the hot component orig- 
inates from both fast and slow winds. Cool filaments that 
enter the inner region appear to be torn by differential rota- 
tion and then shock-heated to the ambient temperature as 
they interact with the hot gas. Since the disc mass is minis- 
cule by AGN standards, the disc is constantly violently af- 
fected by the stars from the other ring. Streams and blobs of 
cool gas are kicked high up above the disc midplane due to 
interactions with winds from stars crossing the disc. These 
blobs rain down on to the disc, and some are brought into 
the sub-arcsecond region of the flow. 



4.2 Accretion on to Sgr A*. 

The 'accretion rate', Mbh, on to Sgr A* is defined in our sim- 
ulations as the rate at which SPH particles enter the sphere 
with radius 7?i n , and is plotted versus time in Fig. 2. Solid 
and dotted lines correspond to time bins of 10 and 200 years, 
respectively. The latter seems to be a reasonable estimate 
for the viscous time-scale of the hot accretion flow at around 
i?in. The average rate we obtain, a few x 10 -6 Mq year -1 , is 
in good agreement with the Bondi estimate (Baganoff et al., 
2003), and is one to two orders of magnitude lower than 
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Figure 2. Accretion rate on to Sgr A* versus time, averaged in 
bins of 10 years (solid) and 200 years (dotted). Both are expected 
to be only rough indicators of the actual accretion rate since much 
of the gas captured may still be unbounded inside R- m and / or has 
to be processed in the inner accretion flow that we do not resolve. 

what previous studies found (cf. §1). However, factors of 
a few variability in the accretion rate are immediately ob- 
vious. This variability can be tracked down to the arrival 
of individual cool gas blobs or filaments in the innermost 
region. Since the density of the hot gas in the vicinity of 
Sgr A* is never zero, the accretion rate never decreases to 
zero. 

It must be stressed that the real Mbh further depends 
on intricate physical details of the inner accretion flow that 
we cannot resolve here. Some of the gas entering R < R[ n is 
unbound and some may become unbound later on as result 
of viscous heating in the flow (Blandford & Begelman, 1999) . 
Therefore, the accretion rate measured in our simulations is 
best understood to yield the outer boundary conditions for 
the inner accretion flow, and as such should be a more phys- 
ically complete estimate of that than the commonly used 
Bondi accretion rate based on Chandra observations of hot 
X-ray emitting gas only. 

The number of particles in the inner 1" (comparable to 
the Bondi radius estimate for this problem) is about 70,000. 
This ensures that we have enough resolution at the inner 
boundary. On average, ~ 4 SPH particles are accreted at 
each time-step. 

4.3 Existence of cold disc in sub-arcsecond region 
of Sgr A*. 

The cold 'disc' found in our simulations is constantly be- 
ing created but it is also likely to be destroyed from time 
to time. Its mass can be estimated as Mdi sc = Mcwt = 
10 M© (M cw /10~ 3 Moyear- 1 ) i 4 , where U = i/10 4 years 
and M cw is the mass outflow rate of cool stellar winds. A 



supernova occurring within the inner 0.5 pc of the Galaxy 
would easily destroy such a disc. The number of young early 
type stars in the Sgr A* cluster is likely to be in the hun- 
dreds, and some have already reached the WR stage (e.g. 
Paumard et al., 2001; Genzel et al., 2003b). Assuming that 
the most recent star formation event occurred a few million 
years ago, one would then estimate the SN rate in Sgr A*'s 
inner star cluster to be at least ~ 10 -4 per year. For com- 
parison, the viscous time at a radius of 1" is 10 6 — 10 7 years 
at best (e.g. fig. 2 in Nayakshin & Cuadra, 2005). Further- 
more, there is much more cold material on scales of sev- 
eral arcseconds and beyond that may be plunging on to 
Sgr A* (Paumard et al., 2004). Arrival of this mass in the 
inner arcsecond could also destroy the disc. 

Therefore it is likely that the cool disc of ~ arcsecond 
scales does not smoothly extend inside the sub-arcsecond 
region of Sgr A*. However, occasionally we observe cool gas 
blobs to directly fall into the capture 'sphere' in the simu- 
lation. It is also conceivable that events started by a super- 
nova shell passage, or by an infall of additional cool mate- 
rial, could also leave a remnant in the form of a cold disc. 
Thermal conduction between the two phases could lead to 
evaporation of a cold accretion disc via the Meyer & Meyer- 
Hofmeister (1994) mechanism, which is unfortunately model 
dependent due to our poor knowledge of the magnetic field 
geometry and strength. 

Nayakshin & Sunyaev (2003) suggested, mainly based 
on the presence of X-ray flares in Sgr A*, that there is a cool 
disc at ~ 0.01-0.1" scales and beyond. The required mass 
of the disc was estimated to be smaller than a Solar mass. 
However, observations of NIR flares (Genzel et al., 2003a) 
put the flares at no more than a few milli- arcseconds away 
from the radio position of Sgr A*, which is somewhat prob- 
lematic for this model. Also, the NIR flare spectra appear to 
strongly favour a synchrotron (i.e. SMBH jet) origin. Fur- 
ther, no eclipses or star 'brightenings' expected when bright 
stars approach the disc (Cuadra et al., 2003) have been ob- 
served so far. Summarising, there is currently no observa- 
tional motivation to favour the existence of such a disc in 
Sgr A*, although it is also difficult to rule out its presence 
(unless the disc extends to the SMBH horizon; Falcke & 
Melia, 1997). 

Concluding, we suggest that there cannot be a smooth 
transition of the larger scale cool disc into the inner sub- 
arcsecond regions of Sgr A*. Nevertheless, the issue of the 
current existence of a cool disc or its periodic appearance 
and disappearance is an open subject for future work. 

4.4 Note on the importance of initial conditions. 

While this work appears to be the most detailed (to date) 
numerical attempt to model the accretion of stellar winds 
on to Sgr A*, observational uncertainties in the stellar or- 
bits and wind mass loss rates and velocities still leave a lot of 
room for uncertainties in the final results. The latest obser- 
vations (R. Genzel & T. Paumard, priv. comm.) reveal that 
the narrow emission line stars might be more equally divided 
between the two discs than what we have assumed here. In 
this case, the disc-like structure would form at a larger scale 
and be probably not as conspicuous. Similarly, if the mass 
loss rate of the 'LBV stars is smaller, the cool disc becomes 
obviously less massive. In addition, mass-loss rates of LBV 
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stars have been observed to vary by more than an order of 
magnitude within a few years (e.g., Leitherer, 1997). This 
effect would bring further variability and uncertainty to the 
results. 

Another important ingredient, still missing in our ap- 
proach, is the inclusion of cooler gas filaments observed to 
be (possibly) infalling on to Sgr A* on several to tens of arc- 
second scales from Sgr A*. These structures, referred to as 
the 'mini-spiral' (Scoville et al., 2003; Paumard et al., 2004), 
would undoubtfully change some of our results. 

Finally, like all previous numerical investigations (Coker 
& Melia, 1997; Rockefeller et al., 2004; Quataert, 2004), we 
have entirely neglected the contribution of the wind pro- 
ducing stars on orbits inside the central 2". While the most 
important wind sources are all located outside of this region 
(Najarro et al., 1997; Paumard et al., 2001), the weaker OB 
stellar winds inside the inner region may still contribute. For 
example, if their wind loss rates are about the inferred accre- 
tion rate, these stars (that appear to be on more randomly 
oriented orbits, see Eisenhauer et al., 2005) could potentially 
destroy the disc in the inner arcsecond. 2 

Future observations of the wind properties and stellar 
orbits near Sgr A* are key to produce an increasingly more 
realistic model of the accretion flow on to Sgr A* . 



5 DISCUSSION AND CONCLUSIONS 

AGN accretion discs are less well understood than discs in X- 
ray binaries on comparable relative scales, because we do not 
have good observational constraints on the origin of gas ac- 
creting on to the SMBH. Sgr A* is becoming the only excep- 
tion to this as observations improve. In this paper we made 
an attempt to realistically simulate the outer ~ 0.1-10" re- 
gion of the gas flow on to Sgr A*. The resulting gas flow 
is far more complex than thought earlier based on studies 
that included non-radiative fast stellar winds from stars ei- 
ther fixed in space or distributed in a spherically-symmetric 
fashion. The presence of cool gas in the sub-arcsecond re- 
gion, as found here, may considerably complicate the in- 
terpretation of observational constraints on the accretion of 
Sgr A*. Although, as discussed above, this does depend 
on still somewhat poorly known details of stellar orbits and 
wind parameters. 

The average accretion rate in our simulation, a few 
xlO~ 6 M0 year -1 , is consistent with the estimates of 
Baganoff et al. (2003) and 1-2 orders of magnitude lower 
than what previous models found (Coker & Melia, 1997; 
Rockefeller et al., 2004; Quataert, 2004). However, this ac- 
cretion rate changes by factors of a few in time-scales shorter 
than hundred years. Then how representative is the cur- 
rent low luminosity state of Sgr A*, if the feeding of the 

2 Note that an extreme case of this was already studied ana- 
lytically assuming a ballistic trajectory approximation by Loeb 
(2004). He showed that if the inner 'S'-stars (e.g., Ghez et al., 
2005; Eisenhauer et al., 2005) produce wind outflow rates as 
strong as <~ 10 — 6 Mq year -1 , then their winds alone could pro- 
vide enough fuel for Sgr A* emission. However, most of the 'S'- 
stars now appear to be of intermediate to later B-type, suggesting 
their winds could be orders of magnitude weaker than assumed 
by Loeb (2004). 



inner region is indeed so turbulent and time- variable as our 
simulation suggests? Afterall, observations of X-ray/7-ray 
echoes from nearby molecular clouds indicate that Sgr A* 
might have been much more luminous some ~ 300 years ago 
(Sunyaev et al., 1993; Koyama et al., 1996; Revnivtsev et al., 
2004). Another aspect of the same issue is that 'accretion' 
of a cool blob in our simulations is not yet a true accretion 
event. If the blob manages to survive in the hot gas and set- 
tles to a disc or a ring at say ~ 10 3 7?s, it may circle Sgr A* 
for a long time without being noticed. Further uncertainty 
in these results is the interaction between the hot and the 
cold gas via thermal conduction. If a cold blob enters the in- 
ner region of the hot flow, and is evaporated there, how will 
this affect the flow there? These and other related questions 
are to be resolved in future work if we want to reach a full 
understanding of the accretion process on to Sgr A* . 
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